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I. INTRODUCTION 


Multiple- scale problems are the norm rather than the excep- 
tion in geophysical flows, and the problem to be described here 
is of this type. In particular, we consider the scenario con- 
sisting of a planetary-scale topography and a baroclinic zonal 
basic flow for which the most unstable perturbations are synop- 
tic-scale. It has been determined, e.g. , by Fredericksen 
(1979), Blackmon (1976), and Blackmon, et.al. (1977), that in 
the Earth's atmosphere, there are regions of preferential 
development of synoptic -scale baroclinic eddies that are depen- 
dent on planetary- scale waves. Fig. 1.1 shows the relationship 
between the Northern Hemisphere oceans and continents , the 
time -averaged SOOmb circulation and bandpass -filtered baro- 
clinic heat- flux data and suggests that the maximum baroclinic 
heat flux occurs just downstream from the trough axes of the 
planetary- scale waves, which are located near the eastern edges 
of the continents . The latitude of maximum heat flux also 
varies, being closer to the equator near the eastern edges of 
the continents than further downstream. 

In the atmosphere, there are two major forcings by which 
the planetary- scale waves may arise: a thermal drive, due to 
land-sea contrasts, and a mechanical drive due to topography. 
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Evidence for the importance of topography in maintaining the 
ultra-long waves is given in Tibaldi and Buzzi (1983); the dis- 
cussion in this study is limited to the effects of topography. 

Early dynamical studies of the effect of topography on a 
flow field were performed for basic states with no zonal flow 
at the surface, and only disturbances of a single wavenumber, 
equal to that of the topography, were considered. These 
studies have resulted in the discovery of baro tropic topo- 
graphic instability by Charney and DeVore (1979), and Hart 
(1979), and extensions to baroclinic situations by, e.g. Char- 
ney and Straus (1980) and Pedlosky (1981). In the study of 
Charney and Straus (1980), a severely truncated two-layer spec- 
tral model with driving and dissipation yields, among the eig- 
ensolutions, a mode which grows in place and other modes which 
are travelling. The height of the topography must be greater 
than a critical value for the stationary growing mode to arise. 
This orographically unstable mode, which draws energy from the 
zonal flow due to its temperature -phase relationship, is pro- 
duced by the form drag and therefore Charney and Straus label 
this instability the "form-drag instability". The travelling 
modes are of the Eady type modified by topography; i.e., they 
are present in the limit of no topography, and are called the 
baroclinic modes . 

Pedlosky (1981) solved the problem of a quasi -resonant flow 
in the presence of topography using weak nonlinearity for both 
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barotropic and baroclinic flows. Stationary Rossby waves are 
resonant with a stationary forcing, which is provided by topog- 
raphy through the form drag (Rambaldi, 1982). Pedlosky (1981) 
considered a flow which is slightly off resonance, in which a 
very small topography is present. He otherwise did not sev- 
erely truncate the model, as in Charney and DeVore (1979), or 
restrict his attention to a strongly anisotropic topography, as 
in Hart (1979), but proceeded with the analysis valid for a 
special ordering relationship for topographic height and dissi- 
pation in terms of a parameter which is small under quasi - 
resonant conditions; namely, the speed that a Rossby wave would 
have if topography were absent. He found that some flows that 
are baroclinically stable in the absence of topography exhibit 
topographic instability. The instability can be obtained for 
arbitrarily small shear, but the width in parameter space of 
the unstable region is inversely proportional to the shear. 

A number of studies have been concerned with the combined 
baroclinic/topographic problem. A recent example is the study 
of Nathan (1985), who has extended the model of Pedlosky (1981) 
by considering a weakly nonlinear problem in which a wavenumber 
n is orographically unstable and wavenumber n+1 is baroclini- 
cally unstable, with comparable growth rates. He demonstrates 
the importance of internal dissipation in determining the non- 
linear evolution of the flow. Attractive features of Nathan's 
analysis are the explicitness of the dynamical relationships 
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obtained and the absence of ad-hoc truncation. Two weaknesses 
of his approach are the requirement for special ordering rela- 
tions between the topography height and other parameters, and 
the lack of scale separation between the orographic and baro- 
clinic wavenumbers. 

i 

The major purpose of this study is to find the linear 
response of the model in the presence of baroclinic instability 
and topography. We will attempt to replicate the parameter set- 
tings for the atmosphere, and more particularly, the annulus. 
We consider the linear baroclinic/orographic instability prob- 
lem in which the coupling of various wavenumbers by the topog- 
raphy can be significant. In contrast, the topographic effect 
on the linear problem in the studies of Pedlosky (1981) and 
Nathan (1985) was not studied explicitly, but brought in at the 
same order as weakly nonlinear effects. The eigenvalue for our 
problem is a complex frequency, and the eigenfunction can be 
approximated by a truncated sum of components of various wave- 
numbers. Linear investigations of this type of flow have been 
done by DeSzoeke (1975, 1983) and by Durney (1977) for an infi- 
nite f -plane or beta plane, with application to the ocean 
mesoscale dynamics. Although the atmosphere allows significant 
meridional wave propagation, we find it useful to consider the 
case in which the baroclinic instability is confined; hence, we 
choose a channel domain. The atmospheric studies cited above 
include a sin Try topography. Our analysis is simplified, in 
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order to bring out the relevant physical mechanisms, if the 
topography is independent of the meridional coordinate . 

The channel model results also are compared with the 
results of the recent laboratory experiments of Li, et al 
(1986), to be referred to hereafter as LKP. Their discussion 
focuses on an experiment run for a choice of thermal Rossby 
number and Taylor number for which wave number 6 dominates the 
flow field when topography is absent. The inclusion of radial- 
ly- independent , wavenumber 2 topography results in the oscilla- 
tion of wavenumbers 4 and 6 at a single frequency. Fig. 1.2 
also shows a wavenumber 2 disturbance whose phase position of 
the trough or ridge oscillates from -15 deg to -45 deg and 
whose amplitude is nearly constant, but exhibiting a frequency 
twice that of the phase oscillation. The theoretical results 
to be presented here show a somewhat smaller mean displacement 
of the trough or ridge and also a somewhat smaller excursion. 
The tendency for strong cyclones downstream of the topographic 
ridge, with high pressure dominating upstream of the topogra- 
phy, was also reported by LKP. The results to be presented 
here are consistent with this finding. The streamfunction 
fields of LKP are describable kinematically as the sum of a 
stationary wavenumber 2 component forced by the topography and 
a travelling disturbance of dominant wavenxamber 4, which con- 
tains also some wavenumber 2 and wavenumber 6 dependence. The 
oscillation of the phase of the waveniimber 2 field requires the 
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Figure 1.2 Experimental time series of amplitude and phase 
from rotating annulus experiments with y- independent topogra- 
phy, after Li, et al. (1986). 
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amplitude of the stationary wave to be larger than the wavenum- 
ber 2 component of the travelling wave. 

The hypothesis of this study is that the flow field of LKP 
is the sum of a stationary wavenumber 2 component and a time- 
dependent wavenumber 4 component. The stationary wavenumber 2 
component could result either from topographic forcing due to 
the presence of a nonzero basic -state zonal velocity at the 
height of the topography, or it could arise spontaneously due 
to topographic instability. The latter scenario would require 
the stationary wave to survive in a finite -amplitude equilibra- 
tion. No unstable stationary eigenmodes are found for the 
relevant region in parameter space; hence the model formulation 
allows for a nonzero zonal velocity in the lower layer, result- 
ing in a stationary forced wave. The time -dependent distur- 
bance observed in the experiments is presumed to be the finite - 
amplitude equilibration of a single eigenmode of a linear baro- 
clinic problem, whose dominant wavenumber is 4. 

In the parameter ranges of relevance to the atmosphere and 
the experiments of LKP, the topography introduces only a 
"slight" modification of the time -dependent flow field; hence, 
this problem can be solved by perturbation methods. The devel- 
opment of the so-called asymptotic model and its results are 
described in chapter 3. A numerical eigenvalue approach, 
described in chapter 4, is valid for flows in which the topog- 
raphy effect is larger. The asymptotic model is found to agree 
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well with the nvimerical model concerning the growth rate and 
the spatial structure of the most unstable modes. 

The discussion focuses on the hypothesized reconstruction 
of the flow fields of LKP in terms of the stationary, forced 
and travelling, free disturbances , and on the comparison of the 
results of the numerical eigenvalue method with the asymptotic 
model results. 



II . FORMULATION AND STATIONARY WAVE SOLUTION 


We used the two- layer quasigeostrophic model to investigate 
the Boussinesq flow over topography in a beta-plane channel. 
The problem is formulated for the case of an arbitrary lower - 
layer basic-state velocity, which induces a forced wave if it 
is nonzero. The assumption of quasigeostrophy implies that the 
Rossby number is small. A second small parameter related to 
the topographic height could be introduced; we take the nondi- 
mensional topography height to be of order Rossby number. 

The quasi -geostrophic equations are developed in the clas- 
sical fashion (Pedlosky, 1979; chapter 6), The scaling is as 
follows : 

(x*,y*) - L(x,y), 
t^ * b/U^ t , 

2 * - H z, 

e* - 0s(z) (1 + Ro Fg 0), (2.1) 

p* - Ps(z) + P(x,y,z,t), 

- Ps ^e P(x,y,z,t), 

/3 - ^0 LV U*, 

where * denotes dimensional variables. The channel width is L, 
and a wavenumber one wavelength is 2nL, The (external) rota- 
tional Froude number and the Rossby number are defined by 
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Fe - fo*LVgH, 

Ro - U*/foL , (2.2) 

and the Coriolis parameter as 

[fo + ^o(y*-V2)]/fo - 1 + Ro ^ (y-1/2). 

Dissipation is present in the form of top and bottom Ekman 
layers , with Ekman numbers given by 

En ~ 1 t H2 n-1,2 (2.3) 

where are the eddy diffusivity coefficients. The 

dissipative coefficients r^j are 

^n - Fn^/V(2 Ro). (2.4) 

We choose our parameters to be of relevance to the atmo- 
sphere, as in Table 2.1. 

Gates (1961) gives numerical values for - g0"^(A0/D) for 
the layer 750 mb to 250 mb, where D— H/2 is the depth of a 
layer. If we average the data from his Table 1 from 1000 mb to 
200 mb, we obtain a value for N of .019 (.013) sec*^ for the 
winter (summer) . The Rossby deformation radius Lq, given by 

Ld - ND/fo, (2.5) 

equals 700 km in winter and 570 km in siammer. Then the Froude 
number F — (L/Lq)^ equals 40 in winter and 60 in summer. In the 


12 


Table 2.1: Parameter values for the planetary scale 


Parameter Atmospheric Value Annulus Value 


L 

4.5 

X 

10« m 

7.5 

cm 

u* 

10 

m 

sec*^ 

.14 

cm sec’^ 

H 

10* 


m 

00 

to 

cm 

D 

5 X 

103 a 

4.2 

cm 

fo 

10-^ 


sec“^ 

4.09 

sec*^ 

S 

10 


-2 

m sec ^ 

10 m 

-2 
sec ^ 

^0 

2 X 

10 

m*^sec"^ 

0 m* 

^ sec*^ 
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annulus, the Brunt-Vaisala frequency for an imposed 2 deg K 
temperature difference and with the thermal expansivity coeffi- 
cient equal to 1.05 x 10’^ (deg K)’^ is found to be N - .5 
sec’^ The resulting deformation radius is about 1 cm, for 
which F - 56. Since wavenumber 4 is observed to be the dom- 
inant wavenumber of the flow with topography, we will consider 
the case of F-35 for which the minimum critical shear is found 
at k=“4 when topography is absent. 

If the Ekman number for the atmosphere is taken to be .001 
(.0001), then V 2 is approximately .65 (.2). In the annulus, 
the Ekman ntimber at the point of interest in parameter space is 
equal to about 2.2 x 10"^, yielding r 2 =- 1.5. The experiments 
of LKP had a free upper surface, in which case r 2 would be 
zero. However, for simplicity, we set the upper and lower dis- 
sipations equal. 

The governing equations, in non-dimensional form, are the 
conservation of quasi-geostrophic potential vorticity equations 
in the two layers (Pedlosky (1979), chapter 6), 

Dt[^^Pl-F(Pi*P2)+;9y] - 

Dc[^^P2‘^(P2‘Pl^‘'’^y5 + J(P2>'?) - -r2V2p2, (2.6) 


In eq. (2.6) p]_ and p 2 , the pressures at each level, acting as 
streamfunctions , are of the leading order in Rossby number; the 
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cotal time derivative is defined as 
^t ^ ^t Pnx^y ■ Pny^x» 

with the appropriate value of n being taken in each layer. 

The (internal) Froude number is defined by 

F - fo^LVN^D^, (2.7) 

where N is the Brunt -Vaisala frequency, given by 

N2 - gAp/(pD), (2.8) 

and the nondimens ional topography height is defined as 

n - h^/(H Ro). (2.9) 

In general, r; is a function of x and y, but in this study, 
we shall consider the case of tf being a function of x only. 
The harmonic analysis of the Earth's topography by Peixoto, et 
al. (1964), points to a dominance of wavenumber 2; however, 
they find significant contributions by odd wavenumbers as well. 
In this study we restrict our attention to topography of wave- 
number 2; i.e., we write 


?7 - 2 hQ cos(2x) ; 


( 2 . 10 ) 


this choice is motivated, in part, in view of the simplifica- 
tion in the modal interactions and also because the experiments 
of LKP showed much simpler flow patterns in this case. 

The value of hQ for the annulus using (2.9) and (2.10) is 
large, on the order of 20. There are several reasons why one 
might expect the effective topography height to be smaller. The 
topographic slope should be scaled to the average radius of 11 
cm rather than the gap width of 7.5 cm. One could also specu- 
late that cold pools of fluid, decoupled from the rest of the 
fluid, may sit in the topographic troughs. (Detailed thermal 
analyses of the experiments of LKP have not been done.) Also, 
the vortex stretching due to the topography is dependent on the 
velocity at the surface, which is smaller than U*. In view of 
these considerations, it is reasonable to set hg/F to 1/7; 
i.e., ho -5 for the annulus , a value we will later show to be 
suggested by the experimental data. 

The boundary conditions are given by 

Pnx ” ^Pnyt) “ 0 at y-0,1 for n-1,2, (2.11) 

where the ( ) denote an average over x. 

The pressure field is written as a sum of a basic state 
and a perturbation. The basic state is a westerly current in 
each layer. We write p]_ and p 2 as 
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Pi " -Uiy + 

P'2 - -U2y + A®^2^ ■*■ A^(^2^' (2.12) 

in which there is present a stationary wave, represented by 
and of scale , and a travelling wave represented by ^ 

whose scale A^ is infinitesimal. 

Substitution of (2.11) into (2.10) yields the equations 

{V2(AS^3^s+^t^^t).p[.(u^.U2)y+AS(^^s.^2®)+^C(^^t.^2^)]+;9y) 

- -riV2(AS^]_s+^t^^t)^ (2.13) 


and 


t a t+ ( aS^2x®+^^«^2x^ ) 3y+ (U2 - aS^ 2y ® 2y ) ^x ] 

+ V2Vx + J(A^<l>2^+^H2^.n) - -r]_V2(A®^l®+A^^l^) , (2.14) 


with 


^nx^"’^nx^*^^nyt^)“^ y - 0,1. 


(2.15) 


The time- independent equations are 


17 


Ui3x 

+^2x®t)3+F(Ui-U2)] - -riV^<i>i^, (2.16) 


and 


U23x l^H<t>2^)-F(4>2^-<f>l^)] + H^J(</>2^,7H2^-^(4>2^-H^)) 

+ ^2x^ ^ ■'■ ^U 2 /A® hQ sin 2x - 4<^2y^^0 

- (2.17) 


with 

^nx®-0 at y - 0,1. (2.18) 

There is no lateral boundary condition for the zonally- averaged 
steady flow; however, the specification of the basic-state 
zonal velocities would imply is constant across the 
domain. Without loss of generality we choose that constant to 
be zero. 

The relative sizes of the terms of equations (2.16) and 
(2.17) are given in Table 2.2. Term 4 of (2.17) is the for- 
cing, which can balance any of a few other terms. If term 5 is 
important, then <f> 2 ^ is a sum of components of even wavenumbers. 
A coupling of zonal harmonics also occurs when the nonlinear 
term 2 is important. We consider the case where terms 3 and 4 
balance. We asstame that we can ignore term 2 of eq. (2.16) and 
terms 2 and 5 of eq. (2.17). This assumption, which we shall 
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Table 2.2: Relative sizes of the terms of equations (2.16) and 
(2.17), after dividing through by F. Numerical values of terms 


in (2.17) 

are given 

for the case * 1, 

U2 - 

.3, 

ho/F - .07, 

r2/F - .03 

for the atmosphere, and the case * 

1, 

U 2 - .3, 

ho/F - .15, 

r2/F - .03 

for the annulus. 




Term 

(2.16) 

(2.17) atmosphere 

annulus 

1. 

Ul 

U2 

.4 


.3 

2. 

A® 

A® 

.2 


.2 

3. 

U 1 -U 2 

U 1 -U 2 

1 


1 

4. 

ri/F 

4U2ho/(FA®) 

1 


1 

5. 


4 ho/F 

.4 


.6 

6. 


r2/F 

.02 


.025 
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call Assiimption A, is strictly valid only when is small. 
Although T 2 /F Is small, r2 is important for determining the 
phase, which is determined by the relationship between dissipa- 
tion and advection by the basic-state velocity. 

If the stationary wave is written as follows: 


“ (ho/F)R®{^n+(y)®’^P(2ix)+bn.(y)exp(-2ix)) , (2.19) 


then we obtain two pairs of equations : 


‘d2 

— + 71+ 
dy2 


^ 1 ± 72 ± ^ 2 ± - 


■d2 

— + 73+ 
dy2 


^2± + 74+ bi+ - -74±. 


with 


bi± “ °2± 


bo+ — 0 


at y - 0,1, where 


( 2 . 20 ) 


( 2 . 21 ) 


( 2 . 22 ) 


71+ - ()3 - FU2 - 4Ui ± 2iri )/(Ui T iri/2) , 
72± - U^F / (Ui iri/2) , 

73+ - (^ - FUi - 4U2 ± 2ir2 )/(U2 T ir2/2) , 
74± - U 2 F / (U 2 =F ir2/2) . 


(2.23) 
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We combine equations (2.20) and (2.21) into a 
fourth-order equation: 


d* d2 

— +(71++73±) — + (71+73+*72±74+) 
dy^ dy2 


b2+ - 0, 


C 


with boundary conditions 


b 2 ± - 


■ d 2 

— + 73+ 
dy^ 


^ 2 ± ">^4± ” ® 


(^ 


at y - 0, 1. The solution is given by 


b2+ “ ^1+ 
b 2 - ” cos 


A++(y-l/2) + A 2 + cos A+.(y-l/2) + A 3 +, 
A.+(y-l/2) + A 2 + cos A..(y-l/2) + A 3 ., 


where 


^++^ - (71++73+)/2 + [(71+-73+)^A + 72+74+]^^^' 
A+.2 - (71++73+)/2 - [ (71+-73+) + 72+74+]^^^* 


A . +2 - (71-+73-)/2 + [ (71-*73.)^/^ + 72-74-]^/^. 

A. .2 - (7i.+ 73_)/2 - ( (7l-*73-)^/^ + 72-74-]^^^. 


single 


.24) 


.25) 


.26) 


.27) 


and where 


21 


^1+ * ("74+ ” ^3+^+- f (A^.^-A^^)] 

^2+ - (A3+A+.^ + 74+)/[<^Os(A+./2) (A+.2.;V^2)]^ 

Ai, - (-74- - A3,A_2)/^eos(A.+/2) 

A 2 - - (A 3 ,A _2 + 74-)/[cos(A../2) (A_ 2 .A,^_ 2 ) j ^ (2.28) 

and 

A 3 + - 74±71±/(71±73±-72±74±) • 

The plots of <^^®/(hQ/F) for F - 35, - 1,5, shown in Fig, 

2 , 1 , for n equal to 1 and 2 , feature an upstream shift of the 
ridge (trough) of <f>^^ from the topographic ridge (trough) 
positions. The phase shift is -10(-15) degrees for the lower 
(upper) layers, a somewhat smaller phase shift than the average 
position of the wavenumber 2 component of LKP, of the wavenxim- 
her 2 component in LKP, 

The terms in equations (2 , 13) - (2 . 15) that are linear in 
yield the following set of equations: 

[at+Ui3x]qi^+ASj(^lS,q^t) + 

[at+U23x]q2’'+A®J(^2®.q2'^) + (U 1 -U 2 ) ] 

+ A®J(^2^*q2^) “ -r2V2^2^. 

where stationary and time dependent potential vorticities 
are defined by 


(2.29) 

(2.30) 
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1 

trough 

i 


ridge 




trough 





Figure 2.1 The order hQ/F stationary flow of wave number 2 
resulting from the interaction of the lower layer basic- state 
zonal flow with topography, for U]^=l, U2=.15, F=35, and 
ri=-r2=1.5. 


i 

! 
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q^s,t ^ [v2^^s , t) ] 

q2^»^ a [V^<^2^ ' ^ ’ (2,31) 

We define the complex streamfunction amplitude by 

<f>^^ - Re{V^^exp(-iwt)}, (2.32) 

for n-1,2. The lower velocity U 2 is taken to be order 1, since 
if U 2 were less than order 1, the amplitude of would not 
be large enough to cause the oscillation in the wavenumber 2 
phase obtained by LKP. Then, the governing equations and 
boundary conditions become 

[ ( - iu>+Uiax) (V2 -F)+riV2+[;3+F(Ui-U2) ] dy^]i>i 
+ (ho/F) [J(^i®,(V 2.F)^1) + J(^i,qi^)] 

+ [-io>+Uiax]F^2 + (ho/F)J(^lS_^2) - 0, 

(2.33) 

and 

[ (-ic^U23x) (V2-F)+r22V2+[;9.F(Ui-U2) ]3 x1t^2 

+ (ho/F) + J(^2><52®>] + J(t^>2>»?) 

+ [-iw+U23x]FV.i + (ho/F)J(^2®>^l> ” 

(2.34) 

with 


^nx"<^nyt>“0 y-0,1. 


(2.35) 
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It is convenient to eliminate in favor of ^2 obtain a 
single equation for ^2> viz., 

{ [ ( -iw+U2ax) (V2 -F)+r2V2+[/9-F(Ui-U2] 

+ (ho/F)[J(^2s.(^^-f)-)-J(q2®. •)]] 

[ ( - iuH-U^ax) (V2 -F)+riV2+[;9+F(Ui-U2) ] 3x (2.36) 

+ (ho/F)[J(fliiS,(V2.F)--J(q2®. •)]] 

+ (-iw+U2ax+(liO/F)J(4^2®.-)] [*i<^Uiax+(ho/F)J(?iiS^.)]F2} 

+ [ ( - iw+U^Sx) (V2 -F)+riV2+[^+F(Ui-U2) ] dy^ 
+(ho/F)[J(fliiS^(V2.F)-)-J(flil®,-)]] (ho/F)4Fsin 2x ^2y “ 0 
with 


^nx-<^nyt>-0 at y-0,1. (2.37) 

The standing wave complicates the analysis by introducing 
y- gradients into the equation for If it were not for the 

boundary conditions, would by y- independent . The plots of 

<j>^ in Fig. 2.1 show the y-gradients to be the largest near the 
boundary. We parameterize <f>^ by considering it to be 
independent of y. We take the amplitude of <t>^ to be .9 times 
the midstream amplitude given by (2.26): 


- -9 [bn+(l/2)exp(2ix) + (b„. (l/2)exp(-2ix) ] . 


(2.38) 
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The choices of the parameters hQ/F and U 2 for subsequent 
analysis are obtained using the experimental data. The lower 
layer velocity scale from the experiments is about .04 cm 
sec“^. The specified ratio U 2 /U]_ is listed in row 1 of Table 
2.3. For these specified ratios, the magnitude of is 

obtained from the stationary wave solution. The terms 
and A^^2^ obtained from the observations, and A^^2^/(^0/^^ ) 
is obtained from the stationary wave solution. The upper layer 
velocity is computed and these quantities are rescaled so that 
the upper layer velocity is one unit. These rescaled quanti- 
ties are listed in rows 3-5. The rescaled observed amplitude 
excursion is shown in row 6. The parameter hg/F is determined 
and listed in row 7. The anticipated wavenumber 2 phase oscil- 
lation, obtained from the theory, is of order (ho/F)A^<^2^» 
is listed in row 8. Note that a comparison of the observed 
versus theoretical values of the phase oscillation; i.e., rows 
6 and 8, would yield agreement for U 2 -.I 5 , for which hg/F is 
equal to .15. 

The excursion in phase, obtained by means of a kinematic 
addition of a standing and time -dependent wave 2, 4 can be com- 
pared with the experimental results shown in Fig. 1.2. If we 
write 

<l>2 * A^ cos (2x+^®) + A^(hg/F) cos[2x-o)^t+^2] » (2.39) 


then we find that 
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Table 2.3: The magnitudes of some observed and theoretical 

quantities as a function of U 2 /Uj_. 


1. ^2/^1 (specified) .1 

2. aS^ 2 s/(^F) (theor.) .25 

3. Rescaled h®<^2s .015 

4. Rescaled A^4>2^ .02 

5. Rescaled aS,^2s/(W^) 

6. Rescaled transient 

wavenumber 2 (obs.) .003 

7. Derived ho/F .2 

8. Hq/F A^<f> 2 ^ (derived) .004 


2 

.3 

.4 

.5 

5 

.8 

1 — [ 
r-t 

1.5 

03 

.05 

.06 

.09 

04 

.06 

.08 

.11 

3 

.8 

1.5 

2.5 

007 

.01 

.015 

CM 

o 

1 

.07 

.04 

.03 

004 

.004 

.003 

.003 
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4> - cos [2x+xl. 


(2.40) 


where \ is given by 


X - -io^Z/2 + tan 


-1 


AS-A*=(ho/F) 
AS-A*^(ho/F) 


tan (wj.t/2) 


(2.41) 


The maximum of tan x occurs for 


tan X 


A t /"t- /-rt\ 


{(AS)2.[(At(ho/F)]2)l/2 


(2.42) 


which is approximately 1/5. This yields a value of x equal to 
approximately 11 degrees, which is on the order of 25 percent 
smaller than the 15 degree excursion illustrated in Fig. 1.2. 



III. ASYMPTOTIC MODEL AND RESULTS 


We consider the case in which a baroclinic flow is modified 
by a very small wavenumber 2 topography; i.e., a planetery 
scale topography. We start with Eqs. (2.23)-(2.25), assuming 
that the zonal wavemamber of the leading order solution, which 
represents a synoptic -scale baroclinic wave, is different from 
0 or 2. The governing equations are 

{ [ ( - iaH-Ui3x) + ] 

+ [(-iw+Ui3x)F+(ho/F)F^2x®3y]l^2 “ 0. 

{ [(-iw+U2ax)(V2-F)+(;3-FU)ax+ r272] 

+(ho/F) 2x]3y)^2 

+ [(-ioH-U2ax)F + (ho/F)F^lx®3y]^l “ 0- 

(3.1) 


with boundary conditions 

^n-<^n>yt-0 (3-2) 

at y - 0, 1. 

We shall use perturbation methods to solve this set of 
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equations for a small topography. We expand and w in terms 
of hg/F; i.e. , 

V-n - + (ho/F)^n^^^ + (V^)" + ... , 

« - w(0) + (ho/F)c.)(^) + (ho/F)2 w(2) + ... . (3.3) 

The streamfunction represents a wavy flow at leading order 
and at order Hq/F. If the wavenumber of the leading order 
problem is 4, there will be a zonally- averaged flow at order 
(hQ/F)^. This component does not affect ^(2) » 
determination of , the application of the solvability 

condition involves a projection onto the baroclinic wave of 
wavenumber 4. 

The leading order problem is the ordinary baroclinic insta- 
bility problem without topography. For the two- layer Phillips 
(1954) model, that problem is described in Pedlosky (1979), 
chapter 7, and is given by Eq. (2. 16) - (2 . 18) with rz-O. One can 
write, in view of the boundary conditions, 

** A^exp(ikx) sin iTry, n-1,2 (3.4) 

which yields a pair of coupled algebraic equations for the 
amplitudes of the two layers. The simultaneous nontrivial 
solution of these requires a determinant of the coefficients of 
and A 2 to vanish, yielding the dispersion relation: 
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- k[U 2 +l/ 2 {c+-ia' ± [c.2-2iac^. - a'2 + 4 a]V 2 )j, 


where 


(3.5) 


c+ - Ug - 2^(K2+F)/[K2(K2+2F) ] , 

c. - {U2K‘‘(K<-4F2)+4/32F2) V 2 /[k 2 (k 2 + 2 F)] , 

R2 - k2 + i2^2_ (36) 

a' - (K2+F)(ri+r2)/[k(K2+2F)] , 

a - {r]_r2K2-ik[r2(y9-UgK2)+r]_(/3-FUg)] )/[k2(K2+2F)] , 

Ug - Ui - U2- 

Since this is a linear problem, the amplitude is arbitrary 
and we set A 2 *l and A]_ 7 , where 7 is given by 


7 - l+KVF + [r2iK2+(^-FUgk]/(Fw<;0)) . (3.7) 


The phase lag T between the lower and upper layers is 
defined by 


7 » I 7 I exp {-i D . (3.8) 


In the case and r]_-r 2 *r, the critical value of Ug; 
i.e., the value for which the waves are marginal, is given by 


Uc - (2r/k) K/(2F - K2)V2, 


(3.9) 
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ucrit,h=0 



0 . 2 . 4 . 6 8 

k 


Figure 3.1 The critical shear, Uj,/2r, as a function of 
zonal wavenumber, without topography, for F=35 and U 2 = .15. 
The minimum critical shear and the short wave cutoff are 
labelled. 
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for which the square root term in eq. (3.4) balances the nega- 
tive imaginary term proportional to r that results from the 
Ekman dissipation. Fig. 3.1 shows a plot of Ug/2r as a func- 
tion of k, which is valid only if 2F is greater than £^ir^ . If 
k is greater than (2F-i^7r^)^/^, no amount of shear will result 
in an instability. A minimvim critical shear of (7(2F) -iw) is 
realized for k*-i?r(y(2F) -iw) . For F - 40(60) and i - 1, the 
minimum critical shear is found at k - 4(5). The positive root 
in Eq. (3.4) represents a perturbation with a westward tilt of 
the troughs and ridges with height, as shown in Fig. 3.2 (for 
t-0, k-5, /3-0, F-35, and rj^ - r 2 ” 1.5), which allows the 

perturbation to extract energy from the basic state. Marginal 
stability occurs when dissipation is large enough to balance 
the energy extracted from the basic state. We select k-4 and 
i-1 to represent the synoptic -scale wave. 

At order hQ/F, we have the governing equations 

-iw(l) [V2 -F]V'i^°^ -io,(l)FiA2^°^ 

+ [ ( - io, ( 0 ) +U]_a^) ( V2 - F) +r [^+F (Ui -U2 ) ] ] ^1 

+ [J(«il®.(V2-F)T6i(0)) + J(|&i<0>,qiS)] 

+ [-iw( 0 )+U]_ax]Fi 62 <^^) + 

(3.10) 


and 


Figure 3.2 The streamfunction for k-5 without topography 
with - 0, F - 35, and r]^ - t 2 “ 1.5, at time t-0, for the two 
layers , with - 1 and U 2 - . 15 . The region shown is the 
half-domain -jr/2<x<7r/2 , 0<y<l. 
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+ [ ( - iw(0)+U2ax) (V2 -F)+r22V2+[/3.F(Ui-U2) ] 3x1^2 
+ [J(?i2®,(V2-F)^2^°^) + 3(^2^°^ •‘12®)] + V> 2 y^°^sin 2x 

+ [-i(j(0)+U2ax]f^l‘'^^ + - 0, 


at y-0 , 1 . 


(3.12) 


The order ho/F governing equations are composed of three 
terms: 1) the leading order operators acting on the perturba- 
tion order ho/F operators acting on which 

is proportional to the order ho/F eigenvalue correction, and 3) 
the order ho/F operators due to the topography and the standing 
wave, which act on the leading order solution. 

In order to find we apply a solvability condition; 

specifically, that any solution to this equation must be ortho- 
gonal to the adjoint of the solution of the leading order prob- 
lem. One can see intuitively, and it is shown formally in 
Appendix A, that is zero, since the replacement of hQ by 

-ho in the equations merely shifts the phase of the topography 
by 180 degrees . 

With - 0, after combining Eqs , (3.10) and (3.11) into 

a single equation, we obtain a forced equation for 
the solution ^2^^^ form 
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^2^^)-g2+(y)exp(i(k+2)x) + g 2 _ (y)exp(i(k-2)x) . (3.13) 


One can solve for g 2 + by substituting 3^ - i(k±2) into (3.10) 
and (3.11), solving in terms of the lower-layer streamfunction, 
and obtaining a pair of fourth-order ordinary differential 
equations in y: 


d‘*g2± 

+ U1+ + ®2±S2± ” “3± (3.14) 

dy^ dy^ 


I 

! with boundary conditions 


•^*621 . 

E2± ~ “4±S2± * (*1)'^“5± “ 0 (3.15) 

dy2 


at y — 0,1, where the coefficients are defined as 

ai+ - -{[(w(0)-U2k)[(k±2)2+F]-hr2i(k±2)2+(k±2)()9-FUs)) 

/(w(0).u2(k±2)^-r2i) 

- {[w(0)-Ui(k±2)H(k±2)2+F]+rii(k±2)2 

+ (k±2)(y3+FUs)}/(w(0)-Ui(k±2)+rii), 
c*2± -[{ [w<0)-Ui(k±2)] [(k±2)2+F]+(k±2)(;3+FUs)-^rl(k±2)2i} 

• { (w<0) ,u2k) [ (k±2)2+F]+(k±2) (^-FUs)+r2i(k±2)2 } 

- {w(0).u2(k±2)Uw(0>-Ui(k±2))F2] 

/[ {o>(0)-Ui(k±2)+rxi) {«^°^+r2i) ] , 
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a3+ - ±2ijr[F{[w(0>-Ui(ld:2)][(ld:2)2+i2«-2+F]+(k±2)(/3+FUs) 

+rj^i[ (k±2)2+i2^2 j )+(lc2+i2^2 

/ { -U(k±2)+r][i] [w^^^+r2 i] ) , 

04 + - -{[w<0)-U2(k±2)]((k±2)2+F]+r2iOd:2)2+(/3-FUs)(k±2)) 

/{w<0)-U2(k±2)+r2i}, 

05 + - T2jri[(4-K2+7F)b2+-F(bi++D] 

/[„(0).u2(k±2)+r2i] . 


(3.16) 


The solutions to these fourth-order ODE's are given by 


S2+ “ ai+cos A++(y-l/2) + a 2 +sin A^..(y-l/2) 
+ a 3 +cos A+. (y-1/2) + a 4 +sin A+.(y-l/2) 
+ ^5+ cos i^ry, 


S2- “ ^ 1 -Cos A^^(y-l/2) + a 2 .sin A..(y-l/2) 
+ a 3 ,cos A..(y-l/2) + a 4 ^sin A..(y-l/2) 
+ a 5 . cos i^ry, 


(3.17) 


where the a^ are constants which must be determined and the 
A's are roots of the characteristic quartic equation, which can 


be found as : 
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A++2 . „i+/2 + [ai+^/h - 02+] V2 

A+.2 - ai+/2 - [ai+2/4 * a2+l^/^ (3.18) 

A. +2 - ai./2 + [ax_2/4 . 

A. .2 - ai./2 - [ai.2/4 . a2.]V2 , 

We can find g]^+ from g2+ by substituting into the lower layer 
governing equation (2.17). 

The functions gj^+ and gjj. are symmetric about y - 1/2, 
which requires that aj^+ - a^^. - 33+ - a3_ - 0 for Jl odd and a2+ 
- a2_ - a4+ - a4_ =» 0 for i even. The specific details of the 

coefficients a^ are given in Appendix B. 

The time -dependent streamfunction field (x,y) is a 

superposition of Fourier components of zonal wavenumber k±2. A 
plot is given in Fig. (3.3) for k-5, i-1, for 0<x<7r; i.e., for 
one "period" for the topography. The functions g2+ have a node 
at y - 1/2, because there. 

It is instructive to write the functions g^+'s as 

gn±(y) “ Gn+(y) exp{i ^„+(y) ) . (3.19) 

It was found that are nearly independent of y, varying 
approximately .OSjt across the y-domain. The upper and lower 
streamfunctions <t>i and 4>2 constructed to order hg/F and are 
given by 
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1^1 - exp(wj^c) { sin iiry cos(kx-w^t-r) 

+ Hq/F G]^+(y)cos( (k+2)x-Wj-C+^]^+] 

+ ho/F G;L.(y)cos[(k-2)x-w^t+^]_.] ) , 

(^2 “ exp(wj^C) {sin iiry cos(kx-Wj-t) 

+ ho/F G 2 +(y)cos[ (k+ 2 )x-Wj-t+^ 2 +] 

+ ho/F G2.(y)cos[ (k-2)x-Wj-t+^2*n . (3.20) 

where ^ Wj. + i 

It is useful to write the order ho flow as a propagating 
wave of wavenumber k that is modulated by a standing wavenumber 
2 pattern. Then, up to order ho/F, the total flow field at a 
fixed latitude can be viewed as an interference pattern between 
the order 1 and order ho flow. The upper streamfunction <f >2 is 
written as 

(f >1 - cos(kx-a)j-t-r) exp[uj^t] {siniiry +ho/F A]^(y) cos [2x*xi(y) ] ) 

+ ho/F sin(kx-w^t - F) expfw^t] Aj^(y) sin[2x*xi(y) ] . 

(3.21) 

where 

Ai = [Gi+2+G^_2+2G3^^Gi.cos({?l++^l.-r)] 1/2^ (3.22) 

and 


^1+-<^1- 

- tan - 

(G 1 +-G 1 .) 

tan 

■ ^i++^i- ■ 

K (3.23) 

2 


[ (G]^^.+G2.) 

2 

J 
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and where F is defined in eq. (3.8). Likewise, the lower layer 
streamfunction is given by 

<i >2 - cos(kx-«j-t ) exp[wj^t] {sinijry +hQ/F A 2 (y) cos [2x-X2(y)]) 

+ ho/F sin(kx-aij.t ) exp[wit] A 2 (y) sin[ 2 x-X 2 <y) ] 

(3.24) 

where 

A2 ■ [G2+2+G2.2+2G2+G2.cos(^i++5i.) ] V2_ (3.25) 


and 


X2* 


^2+‘^2- 

f (G 2 +-G 2 .) 

- tan ■< tan 

82++Q2- 

2 

l (G24.-K32-) 

2 


(3.26) 


At longitude 2x - Xn> order hg flow in layer n adds 

constructively to the order 1 flow, and at longitude 2x - 
is found a destructive interference. Since the order hg 
streamfunction is antisymmetric for i - 1, the interference 

pattern for the south half- domain 0<y<l/2 is shifted by jt/ 2 
from the north half- domain l/2<y^l. The functions A^ and Xn 
are plotted in Fig. 3.4. On each side of the mid- channel lati- 
tude, the function Xn nearly independent of y. The lower 
layer amplitude A 2 is a maximum at y-1/4 and y-3/4, and the 
upper layer amplitude maximum is slightly further from the mid- 


channel . 
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Figure 3.4 Amplitude (An) and phase (^n) constructive 
interference between order 1 and order Hq/F solutions, for 
layers n=l and n=2 , as a function of y - 1/2 . For 0<y<l/2 , we 
note An(l/2-y) = ANd/2+y) and Xnd/2-y) = Xn(V2+y) ±^r. 
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The total time -dependent flow field, shown in Fig. 3,5 for 
times * 0, tt/2 , tt, and 37t/2 , is a procession of highs and 
lows whose centers meander as they move east, on a track whose 
position is constant in time. For the lower layer, the furthest 
south excursion is found at about x*7t/3 and the furthest north 
excursion is found at about x«-7t/6. The amplitude is smaller 
for the upper layer, and the furthest north excursion is 
slightly upstream of the topography ridge, and the furthest 
south excursion near the trough. 

One might interpret the avoidance of the low latitudes near 
X =- -tt/S as the presence of a "blocking ridge" , immediately 
upstream of the topography ridge. This feature of the results 
is similar to those of LKP. However, LKP found also a strong 
longitudinal dependence of the strength of the highs and lows; 
only a weak variation in strength was found for the time- 
dependent solution alone. The addition of the stationary wave 
yields a stronger longitude dependence in the amplitude of the 
wave . 

In the atmosphere, Fig. 1.1b, after Blackmon et al. (1977), 
suggests that the main storm track in winter is further north 
upstream of the continents than it is downstream. Although the 
relationship between the path of the cyclones and anticyclones, 
the eddy heat flux, and the time-mean circulation was noted in 
the Introduction, it is interesting that the storm track is 
very evident in Figs. 3.5 and 4.3, while there is very little 
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Figure 3.5 The streamfunction <f>^ for layers n=l,2, 

computed for times = 0, tt/2, tt, and 37t/2, using the 

asymptotic model. Only the half of the domain -7r/2<x<7r/2 is 
shown. Topographic ridge and trough positions are shown by 
arrows . 
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longitude dependence of the intensity of the waves, and the 
most unstable mode, divided by the exponential growth factor, 
has a zero time mean; i.e. the time-averaged circulation is 
simply the zonal basic state. We have demonstrated the possi- 
bility that storms follow a meandering path without the pres- 
ence of a wave in the time -mean circulation. 

At order (hQ/F)^ we find the lowest order eigenvalue correc- 
tion, which is given for the case of U 2 - 0 by 

y-1 

o,(2)i«i/+4^2^[^(0) -u^+ir^] [cos iTry (g.y-g+y)]| 

y-0 

(3.27) 

where I and I' are integrals given in Appendix A. The results 
for arbitrary U 2 are also given in Appendix A. These integrals 
are obtained by writing the order (hQ/F)^ equations and 
applying the solvability condition, for obtaining an expression 
for The streamfunction that order has 

components of wavenumber k, k+4, and k-4; only the component of 
waventunber k contributes. 

The marginal stability curve was computed from this small 
topography perturbation analysis by choosing F and finding the 
critical dissipation using a Newton's iterative method. A 
quantity AU^ which is the percent difference between for 

topography height hg - 5 and for no topography is plotted 





Figure 3 . 6 The percent change in the critical shear as 
a function of dominant wavenumber k due to the presence of 
topography, with F=35, ^9 = 0, r2^=^r2==1.5, and hQ = 5. The 
solid curve is computed using the asjrmptotic model, while the 
circles are obtained from the numerical. No data is obtained 
at x«2 since the as 3 rmptotic model results are invalid there. 
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as a function of k in Fig. 3.6, for a fixed F - 35. A nega- 
tive AU^ implies destabilization. It can be seen in Fig. 3.6 
that the topography destabilizes the flow near the short-wave 
cutoff for each given wavenumber. 

The time -dependent equilibrated amplitude was chosen to be 
equal to .025, the magnitude the rescaled term from table 
2.3 would be when U 2 /U^ is equal to .15. 

The sum of the time -dependent field and the stationary wave 
field, shown in Figure 3.7, reveals a prominent standing wave 
of wavenumber 2, with a less prominent wavenumber 4 pattern 
superimposed on it. The combination of the topography and the 
parameterized standing wave acts as an effective y- independent 
topography in both layers. Physically the meandering of the 
storm tracks is a result of the antisymmetry of the vortex 
stretching about the middle of the channel due to the interac- 
tion of the effective topography with the u- velocities associ- 
ated with the leading-order disturbance. The antisymmetry is 
masked in the total wave field due to the presence of the 
larger stationary wave. 

A tendency toward high pressure immediately upstream of 
the topographic ridge and low pressure downstream is evident 
after the addition of the stationary wave component. The phase 
of this feature relative to the topography agrees qualitatively 
with the similar features of the atmosphere and the laboratory 


experiments . 
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The antisymmetry of the wavenumber 2 time -dependent solu- 
tion about the middle of the channel yields a difficulty in the 
interpretation of the phase oscillation of the wavenumber 2 
component. According to the theory presented here, the phase 
oscillation should be zero at the middle radius and a maximum 
at y-.25 and y*.75. The observed magnitude of the phase oscil- 
lation is nearly independent of the radius. Some of the sim- 
plifications in the theory, for example the lack of a station- 
ary wavenumber 4 component and the parameterization of the sta- 
tionary wave solution may be responsible for this discrepancy. 


IV. NUMERICAL APPROACH AND RESULTS 


4.1 Formulation: 

More general results, i.e., results which are still 
obtained from a truncated model but one in which we relax some 
of the smallness restriction on hQ, can be found by numerical 
methods. The numerical approach in this study was used only 
for the time -dependent problem. The parameterized stationary 
wave solution, forced by the topography in the presence of U 2 , 
defines a coupling of the various Fourier components in a man- 
ner similar to the topography term. In contrast, many other 
studies, for example, Pedlosky (1981), consider only a zero 
lower layer basic state velocity. We use Assumption A 
described in Chapter 2; namely, assuming the stationary wave 
parameterization to be valid for arbitrary hQ/F. In this study 
the main purpose of the numerical approach is for comparison 
with and validation of the asymptotic time -dependent model 
results presented in Chapter 3. 

We define as 

CO CO 

- X cos irry + ^ exp(ikx) 

i-1 k— CO (4,1) 
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and we recall that 

7 ^ - - 4 hQ sin 2x. (^-2) 

We substitute the above expressions into eqs (2.20) and get 
for the zonally- averaged flow the following equations, 


w [ (i2;r2+F)ao^l - Fao^2^ aQ^x 

+2;ri(ho/F){bi+[-(22+i27r2)a.2il+Fa2^2] 

-hX- [ (22+i2,r2)a2^X+Fa2^2] 

+ [ (4+F)bx+-Fb2+]a_2^X [ (^+F)bx. -Fb2. ] a2 ^x3 " 


w [ (i2fl-2+F)ao^2 ■ ^^Oil] ir2i^7r2 aQ^2 

+27ri(hQ/F) {b2+[ (-2^+i^7r2)a.2^2‘^^^2il] 

*^2- f (2^+-827r2)a2^2'*'^^2il^ 

+ [ (4+F)b2+* Fbx+] a. 2.22 [ (^■*"F)b2. "Fbx. ]a2^2] ” ® 

+ 2i;rho[a2.g2'^-2i2] ” 


and for each zonal nonzero wavenumber k, we obtain the follow 
ing equations: 
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[(w-Uik)(K2+F)+k()3+FUs)+irK2]aij^l - (w-Uik)Fak^2 
+2ffiho/F{ , 2 f ^(k- 2) (k- 2) il 1 

+[-( 4 +F)bi++Fb 2 +]a(k. 2 )ill 1 
-5k, -2t^l-t'(^i+^'^^)^(k+2)i2+^^(k-2)il] 
+[-( 4 +F)bi.+Fb 2 .]a(k+ 2 )il ] 1 

imax 

+ 2 iirhoI [Pmif 
2-1 

(l- 5 k, 2 > [bi+[-(K 2 +F)a(^. 2 )m 2 +^^(k- 2 )ml 3 
+ [ - ( 4 +F)bi++Fb 2 +]a(ij. 2 )nil j 

-<l-^k, - 2 ) [l 5 i.[-(K 2 +F)a(k+ 2 )m 2 +^^(k+ 2 )mll 
+ [-( 4 +F)bi.+Fb 2 .]a(i^+ 2 )min ” 

and 


[(a)-U2k)(K2+F)-Hk(/3-FUg)+irK2]ak^2 ’ (‘^-U2k)Fai^.^l - 0, 
+ 25 r 2 ho/F{ 5 i^ 2 t^ 2 +t*(^ 2 -^+^>®(k- 2 ) 21 +^^(k- 2 )i 2 l 
+ [-( 4 +F)b 2 ++F(bi++l)]a(;^. 2 )i 2 l 1 

-5k,.2[t>2-[*(*^2+^+^)®(k+2)21+^^(k-2)22] 

+[ - ( 4 +F)b 2 .+F(bi.+l) ]a(k+ 2 ) 22 l J 

imax 

+2i7rhoX[pm^{ 

2-1 

(l- 5 k, 2 > tb 2 +[-(K 2 +F)a(k. 2 )ml+F^(k- 2 )m 2 l 
+[ - (4+F)b2++F(bi++l) ] a(k.2)m2l 

* ( 1 -^k, - 2 ) tb 2 - [ - (K 2 +F)a(k+ 2 )ml+^®(k+ 2 )m 2 l 

+ [-(4+F)b2.+F(bi.+l)]a(k+2)m2n “0 (^-3) 
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where is defined by 


cos iwy - X Pmi 
m-1 


i.e. , 




' 2A[l/(m+i)+l/(m-i)] 

0 

V 


if m±i odd, 
if ra±i even. 


k' * ^ k-k' ; otherwise 


and 


Kjjj+2 , (k±2)2+m27r2. 


(4.4) 


(4.5) 

(4.6) 


The matrix eigenvalue approach is taken, in which the state 
of the system is represented by Q and the dynamics of the 
system are represented in the form 


Qt + A Q - 0. (4.7) 


If we write Q - Re { q exp(-iwt)}, then eq. (4.7) gives 


(A - i cj I) q - 0. (4.8) 


I 


Define potential vorticities 
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Qkin ” ^ ^kiO-n)* 

for n-1,2. Then we can express the streamf unction amplitudes 

as 

^kil “ *[FQkil+<^^'^^^Qkill/%- 

^ki2 - -[FQki2+(K"+F)Qki2]/KF. 

where Kp-K^ (R2+2F) . Equations (4.3) can be rewritten in terms 
of the Qij-jgi and Qkjj2- 

irii^^2[p Qq^ 2 + Qoiil/Kp 

+ (ho/F)2ffi{ [Q_2il^l+*Q2il^l- ] 

+ [(4+F)bi+-Fb2+][FQ.2i2+(K^+F)Q.2ii]/[22(22+2F)] 

+ [(4+F)bi.-Fb2.][FQ2^2''’(^^+F)Q2ii]/[22(22+2F)] ) 

- " QOil ” O' 

ir2i^»r*(F + (i^w^+F) Qofll/Kp 

+ 2-eirho/F{ [Q_2i2'^2+**^2i2^2-l 

+ t (4+F)b2+-Fbi++F] [ FQ.2il+(K2+F)Q.2i2 1 / 1 2^ (22+2F) ] 


[(4+F)b2.-Fbi.+F][FQ2J!i+(K2+F)Q2i2]/t22(22+2F)]) 

-w Qo^2 “ 0. 
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kUlQkil- [k(^+FUg+iriK2 ] ( FQkj,2+(k^+F)Qkj;i ] /Kp 

imax 

+ 2ijr(ho/F)(5]{Pi„^(l-5k,2) I Q(k+2)ml^l- 

m -1 

- [ (4+F)bi+-Fb2+] [FQ(k-2)m2+(K^+F)Q(k-2)ml] APm- 1 

- Pm2 (l-^k,-2> C Q(k-2)ml^l+ 

- ((4+F)bi.-Fb2.] (FQ(k+2)m2+(K2+F)Q(k+2)ml3/%m+l 

+ ^k,2 [QOil^l+ ■ K^+F)bi+-Fb2+] [FQo^2+(K^+^>QOil]/KFm-] 

■ ^k,-2 tQO-ei^l- ■ [ (^+F)bi.-Fb2.] [FQo^2"'’(^^‘''^)QOil]/KFm+] 

-w - 0, 

and 


kU2Qki2- [k(y9-FUg)+ir2K2] [FQ^j;i+(K2+F)Qkjj2]/KF 

+ 2ix(ho/F)(X{Pmi(l-^k,2) t Q(k+2)m2^2- 
m-l 

• [ (4+F)b2+-F(bi++l) ] [FQ(k.2)ml+(K2+F)Q(k.2)m2]/KFni- ] 

- PmJl (l-^k,-2> [ Q(k-2)m2b2+ 

- [(4+F)b2.-F(bi.+l)] (FQ(k+2)ml+(K2+F)Q(k+2)m2l/KFm+] 
+ ^k,2 [Q0i2'^2+ -K^+F)b2+*F(bi++l)] [FQoii+(K2+F)Qo^2]/KFm-l 

* ^k,-2 [Q0i2^1- * [ ] [FQo^l+(K2+F)Qo^2]/KFm+3 

-w Qy._£2 - 0. 


(4.11) 
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where 

ip « i2^2(_£2^2+p) ^ 

i^FDd:*K^"(Kni±'+2F), 

K2^2.22+i2^2^ 

KF2i-*K2i2 (K2 jj2+2F) . (4 . 12) 

Each equation in (4.11), representing one row of the matrix 
in eq. (4.8), is used to compute the matrix coefficients, and 
the eigenvalues and eigenvectors associated with the matrix are 
found with standard eigenvalue solving mathematical subrou- 
tines . 

The coupling in (4.11) is between the zonal wavenumber k 
and k±2k' . So the equations for odd k are decoupled from even 

k. For odd k, there is no coupling between a wavenumber pair 
(k,i) and (k±2,i), since - 0. For even k this is also 
true, except that (-2,i), (0,i) and (2,i) are all coupled 
because did not enter into these equations. Then, eqs. 
(4.11) can be written for four collections of wavenumbers, each 
consisting of wavenumber pairs that are coupled, and which are 
totally decoupled from wavenumbers present in other collec- 
tions . We shall refer to these collections as set numbers 

l, . . . ,4. By convention, sets 1 and 2 will have even zonal 
wavenumbers and sets 3 and 4 have odd zonal wavenumbers. Set 
number 1 has (k,i) - (0,1) as a member, and set number 3 has 
(k,i) - (1,1) as a member. As an example, if the system of 
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equations is truncated at kmax - 6 and imax - 4; i.e., -6^k<6 
and 1 <1 < 4 , we will define set numbers 1,..,4 as follows: 

Set #1: 

(k.i) € {(-6.1). (-6. 3), (-4. 2). (-4. 4), (-2.1), (-2, 3), (0,1), 

(0,3), (2,1), (2, 3), (4, 2), (4, 4), (6,1), (6, 3)) 

Set #2: 

(k.i) € {(-6. 2), (-6, 4), (-4,1), (-4, 3), (-2, 2), (-2, 4), (0,2), 
(0.4), (2, 2). (2, 4), (4,1), (4, 3), (6, 2), (6, 4)) 

Set #3: 

(k.i) € {(-5, 2), (-5, 4), (-3,1), (-3, 3), (-1.2), (-1,4), (1,1), 

(1.3) . (3. 2). (3, 4), (5,1), (5. 3)) 

Set #4: 

(k.i) € ((-5,1), (-5, 3), (-3, 2). (-3, 4), (-1.1), (-1,3), (1.2), 

(1.4) , (3,1), (3. 3), (5. 2). (5, 4)). 


4 ■ 2 Growth Rate : 

The growth rate is decreased only slightly by the inclusion 
of topography for the parameter ranges of interest. For F-35, 
^1“^2 " 2, the growth rate is decreased on the order of .04 
day'^ as the topography height hQ is increased from 0 to 10. 
We can define a function u(hQ;w^) to be the shear U that yields 
a growth rate of for a topography height hQ. We find that u 
is approximately given by 
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Figure 4.1 The transition from wavenumber 4 to wavenumber 
5 as the dominant zonal wavenumber of the fas test -growing eig- 
enmode (solid line) , and the marginal curve (dashed line) , as 
functions of the shear Ug and topography height hQ, for F = 35, 
U 2 = .15, ^ = 0, and r^^ « r 2 * 1.5. The dominant meridional 
wavenumber is i = 1 . 
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where a is approximately .002. 

Of greater importance are structural changes introduced by 
the topography. Near the minimum critical shear, the as 3 mip- 
totic model showed long waves to be stabilized more than short 
waves (Fig. 3.6), and hence, for some parameter values the 
topography can give rise to an increased dominant wavenumber of 
the fastest growing mode. An increase in the zonal wavenximber 
as ho increases is suggested by Fig. 4.1 for 0<hQ<7 ; and for 
ho>7, the wavenumber decreases, but at this parameter setting, 
the wave is not unstable. 


4.3. Streamfunction: 

The relative amplitudes associated with the various wave- 
numbers for the most unstable eigenmode when F«35, ho *• 5, and 
^1*^2 “2, is shown in Fig. 4.2. The dominant wavenumber is 
(k,i) - (5,1). The sidebands picked up by the asymptotic model 
are of wavenumbers 3 and 7. A wavenumber 1 effect that is not 
part of the asymptotic model solution is also present. In the 
lower layer, this wavenumber 1 component is larger than the 
wavenumber 3 component . 

Time series of the streamfunction field for the above 
choice of parameters are shown in Fig. 4.3, where only half the 
zonal extent of the domain is shown. The streamfunction field 
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amplitude 



Figure 4.2 The relative amplitudes of various Fourier 
components making up the most unstable eigenmode, for F=35 , 
U2*=.15, ^^0, hQ = 5, and ^ 1.5. 



upper 


lower 


Figure 4.3 The streamfunction for layers n=l,2, 

computed for times = 0, tt/2, tt, and 37 t/ 2, using the 

numerical model, for U2=.15, F=35, hQ = 5 , /3==0 , and "" 

1.5. Only the half of the domain -7r/2<x<7r/2 is shown. 
Topographic ridge and trough positions are shown at the top. 
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agrees very well with that obtained from the asymptotic model 
in Fig. 3.5. 


4.4 Effects of Truncation: 

First we consider the inviscid problem in which the Fourier 
expansion in eq. (4.1) is truncated to the zonal flow and the 
topography wavenumber. This is the simplest scenario in which 
we could address the question of whether a stationary wave 
might arise form the case of U 2 equal to zero. A nontrivial 
solution to eqs (4.3) can be found by setting the determinant 
of the matrix of coefficients of the equal to zero. The 

result is the following biquadratic equation: 

[w^+b]2 - c^w^ = 27 [u )2 + b+ac], (4,14) 

where 

2 (/9 + FU) 2 (/9-UK2) 

a = -2 U + = , (4.15) 

K2 + F K2 + F 

4 (/9 - FU) (y9 - U K^) 

b = , (4.16) 

K2 ( K2 + 2 F) 

4 (K2 + F) 

c = - 2 U + , (4.17) 

K2 (K2 + 2F) 


and 
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Figure 4.4 The right-hand side and the left-hand side of 
the dispersion relation, eq. (4.14). The curves shown are for 
F*0 (solid), F-7 (dashed), F-14 (dotted), and F—IOOO (dot- 
dashed) . 
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4 ho^ ( ^ 2^2 + F) (K2 + F) 

7 = . (4.18) 

(i2;r2 + 2F) R2 ( k2 + 2F) 

Although w2 is in general complex, the case of real 
gives us some insight into the nature of the solution. The 
left-hand side of eq. (4,3) is a parabola centered at 
(c 2 / 2 -b, c 2 (b-c2/4) ) and the the right-hand side is a line of 
slope 7 passing through the point (-b-ac,0). A solution is 
obtained whenever the two sides of eq. (4.3), shown in Fig. 
4.4, are equal. If the line and parabola intersect, then there 
is a solution for real . If o>^ is positive, then we have a 
neutral eigensolution. If is negative, the solution will 
be pure imaginary; i.e., the perturbation is growing without an 
oscillatory behavior. If the line and parabola do not inter- 
sect, then the perturbation has a complex growth rate; i.e., it 
has a growing and an oscillatory behavior. 

In Fig. 4.4, we see that for F > 7, without topography 
is complex. The solution is a travelling baroclinic wave with 
k ±2. If F > 21, then as 7 increases, implying increasing hQ , 
there are two critical values of 7 , namely 7 ^]^ < 7^2 such that 
for < 7 < 7^2 there are four neutral solutions, and for 7 

> 7^2 there are two neutral solutions and two non- oscillatory 
growing solutions. For 7 < F < 21, 7^2 does not exist. For F 
less than 7, topography destabilizes neutral modes. The 
remaining discussion is for F > 21. 
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upper 


lower 
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The streamfunction pattern for the case of 7 < 7d» shown 
in Fig. 4.5, is a pattern of highs and lows travelling on a 
meandering path, as was the case of a baroclinic wave modified 
by topography presented in Figs. 3.5 and 4.3. In addition, 
Fig. 4.5 reveals an index cycle vacillation, due to the pres- 
ence of a zonal flow in the perturbation. At time w^t equal to 
0 (tt) , the flow has a high (low) zonal index. 

The case of 7 > 7^2 » shown in Fig. 4.6, is a stationary 
disturbance, with a ridge centered over (shortly upstream from) 
the topographic ridge in the lower (upper) layer. 

The above discussion concerns a severe truncation. Except 
for a restricted region in parameter space, viz. for relatively 
small F, the 2,0, -2 truncation does not yield the most unstable 
mode. Hence this mode may not be physically relevant, unless 
it dominates in a finite amplitude problem, as some less 
unstable modes do in Hart (1981) and Nathan (1985). The numer- 
ical eigenvalue problem was run for the viscous case to deter- 
mine if this mode arises at the parameter setting we have con- 
sidered. The critical topography height at which the topo- 
graphic instability for F-35 is equal to approximately 6 for 
the inviscid case, and increases for viscous cases. Hence it 
is more reasonable to assume that the stationary wave arises 
due to a nonzero lower layer basic state velocity, and not due 
to a topographic instability with a zero basic state lower 
velocity. 



Figure 4.6 As in Fig. 4.5, except for hQ = 6 , and only at 
time t = 0 . 
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Table 4.1 : Growth 

rates as 

a function 

of wavenumber 

selected topography heights, 

for two truncation levels 

Truncation 1: 

Kax - 12 

> -^max ” 6: 


ho: 

0 

5 

6 

Dominant k : 1 

.427 

.216 

.172 

2 

.82 

0.0 

.314(wr”0) 

3 

1.14 

.952 

.896 

4 

1.36 

1.22 

1.18 

5 

1.45 

1.36 

1.32 

6 

1.37 

1.33 

1.31 

7 

1.03 

1.05 

1.05 

8 

0.0 

.29 

.37 

Truncation 2: kmax - 


ho: 

0 

5 

6 


0.0 .334(wj.-0) 


some 


Dominant k : 2 


.82 
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In table 4.1 is shown a plot of growth rate as a function 
of wavenumber for the inviscid dase with U 2 - 0, for hg - 0, 5, 
and 6 , for ^ hQ - 0 , 5 , and 6 . 

Stabilization due to the topography occurred for all wavenum- 
bers, except near the short-wave cutoff and also with the nota- 
ble exception of wavenumber 2, which exhibits the topographic 
instability for hQ ^ 6. Hence the severely truncated scenario 
described above is present when the severe truncation is 
relaxed. 

The effect of truncation on the most unstable eigenvalue cj 
for F - 35 and r 2 ^-r 2 - 1.5 were determined for "" ^max 

10. Results for hQ - 5 are shown in Table 4.2. For fixed 

kjjjax-®» it appears that ^max * 6 is sufficient meridional 

resolution to determine the eigenvalue accurately. Also “ 

8 appears to be sufficient zonal resolution. The appropriate 
value for k^ax» ^ lesser extent >?jnax increase with 

increasing F, since the dominant wavenumber is a function of F 
as given in Chapter 3. The number of wavenumbers required in 
these results is not nearly as great as in Niehaus (1980) , who 
studied wave stability. One might suspect that the coupling of 
zonal harmonics in wave -stability problems would be stronger 
than in this problem, where imposed waves do not occur at lead- 
ing order. 

The change in growth rate, Ac4;^, and frequency, 
introduced by topography were computed for the numerical model 
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Table 4.2: The nondimens ional growth rate as a function of 
the modal truncation, for ho - 5, F - 35, 0 - 0 , U 2 - .15, and 
ri - r2 - 1.5. 


^nax 

■^max 

Growth 

4 

2 

.276 

4 

>4 

.273 

6 

2,4.6 

.303 

6 

>8 

.297 

8 

2 

.273 

8 

>4 

.270 

10 

2 

.273 

10 

2:4 

.270 
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and compared with the asymptotic model. The percent discre- 


pancy was computed using the following definition: 


[(Awi>n - (Awi>a] 

^ni ” > 

[ (Awi)a^+(Awj.)a^] 


[(Awr>n - (Awr)aJ 

^nr ” : ’ ^^0% , 

[(A«i)a2+(Awr)a^]^^^ 


where 


(Awj.)a ” “rl - “rl 

hQ no topo . 


(Awi)a - «i| - Wil 

hQ no topo . 


for the asymptotic model, and 


(Awr)n “ “rl ' “rl 

hQ no topo . 


(Ao>i)n - wil - Wil 

hQ no topo. 


for the numerical model with - 10 and ““ these 

results are plotted in Fig 4.7. Most of these results are 
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percent discrep 



hzero 


Figure 4.7 The percent discrepancy € 5 j- (solid) and 
(dashed) , between the topographic effects on the frequency and 
growth, respectively, for /S=0, F=35, and = ^2 = 1.5, for 

various values of hQ, where l%ax “ 
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almost identical to those obtained with l%ax * residual 
in the discrepancy as k^iax -^max infinity measures the 
error inherent in the small hQ approximation used in the asymp- 
totic model. 

For the point in parameter space given in Fig. 4.7, the two 
models disagree much more in the growth rates than in the 
frequency change Aw^. Also we find the absolute value of Aw^ 
to be smaller than the absolute value Acj^. This is not true 
when one looks at some other points in parameter space; e.g., 
F-20, shown in Table 4.3. In all the cases run, for hQ less 
than or equal to 5 , we find the frequency to be increased by 
the presence of the topography. 

Since from the asymptotic model does not project onto 
the wavenumber of the leading order disturbance, one would 
expect the models to agree to order (hg/F)'*. The errors appear 
to be of order (hg/F)^, due primarily to the fact that coupling 
terms due to the topography in the matrix in equation (4.7) 
include a factor of 2n2: the factor of 2 from the x-derivative 
of the topography height, and nJl from the y-derivative of the 
wave streamf unction. The largest discrepancies occur when F is 
large, which means hQ is large when hQ/F is fixed. The ratios 
of some terms that do not depend on F; for example the compari- 
son between the advection and topography terms in the original 
formulation, might become important. 
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Table 4.3: The percent discrepancies and for various 
values of F/r 2 and F, with hg/F equal to .1, U]_ - 1, U 2 - .15. 


F 

F/V2 



20 

20 

6.75 

4.68 

40 

20 

3.1 

11.1 

60 

20 

.05 

15.4 

20 

200 

8.92 

6.09 

40 

200 

3.57 

17.08 

60 

200 

3.89 

26.67 


4 . 5 Phase Speed : 

The dominant wavenxamber, i.e., the wavenumber k of the 
leading order solution, can be used to define a phase speed; 
i.e. , 

Cj. — Wj./k 

The phase velocity was found to be affected only slightly by 
the topography, increasing on the order of a few per cent as 
the topography is varied from 0 to 5 . 


V. CONCLUSIONS 


The observations of some features of the flow field in the 
atmosphere and in the annulus experiments of LKP have been com- 
pared with results of a linear theory. Notably, the observa- 
tions indicate a tendency for a ridge to exist upstream of the 
topography due to the presence of a stationary wave forced by 
the interaction of the lower layer zonal flow and the topogra- 
phy. 

The transient disturbance involves a coupling of various 
zonal wavenumbers by the topography and the parameterized sta- 
tionary waves. When parameters are chosen appropriately for 
the atmosphere and the annulus, the fastest-growing eigenmode 
is a baroclinic wave, modulated by the topography, causing the 
eastward- travelling cyclones and anticyclones to meander along 
a sinusoidal path that is fixed in time. The cyclone and anti- 
cyclone centers reach their minimum latitude at about one- third 
of a topographic wavelength downstream from the ridge, and the 
maximum latitude one-sixth wavelength upstream. The phase of 
this path relative to the topography is in qualitative agree- 
ment with the observational results of Blackmon, et al. (1977) 
shown in Figure (1.1b). 

The form of the transient disturbance is essentially the 
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same for the cases of zero and nonzero lower layer basic state 
velocity. Hence the meandering wave is present for the case of 
no lower- layer velocity, but in this case the forced stationary 
wave is absent. 

The transient disturbance alone cannot capture much longi- 
tude dependence of the strength of cyclones and anticyclones. 
However, the sum of the stationary and transient disturbances 
yields a flow pattern characterized by a tendency toward high 
pressure shortly upstream of the topographic ridge and low 
pressure shortly upstream of the topographic trough, which is 
in qualitative agreement with the observed results. 

The eigenvalues obtained in this study show the topography 
to stabilize baroclinic waves except near the short-wave cutoff 
and to shift the maximiam instability to shorter wavelengths for 
small topography. The eigenvalue changes are small, of order 
(h0/F)2. 

The asymptotic and numerical model results agree to order 
(hQ /F)^. Hence there is good agreement for moderately large 
hQ, since F is large enough to keep hQ/F small. Physically, 
large F implies that the length scale of the most unstable wave 
is much shorter that the wavenumber 2 topography; hence the 
magnitude of topographically- induced vertical motions on a bar- 
oclinic wave is smaller than for an equivalent problem with 
small F. One might suspect that the agreement should be of 
order (hQ/F)^, based on the as 3 nnptotic model formulation. Some 


76 


further check on the magnitudes of the various terms might be 
required to verify this. 

If one truncates the Fourier expansion of the perturbation 
in the numerical model to include only the zonally- averaged 
flow and the zonal wavenumber of the topography, the inviscid 
dispersion relation reduces to a biquadratic which yields a 
topographically-modified baroclinic instability for small 
topography and a topographic instability, in which waves grow 
in situ, for large topography. For intermediate values of hQ, 
the waves are neutrally stable. The eigenmodes obtained from 
this truncation were not the most unstable for the values of F 
relevant to the planetary-scale atmosphere. It remains to be 
seen if they might be important in a nonlinear model. 

This study can be extended in several ways. The stationary 
wave problem was arbitrarily simplified, in order to facilitate 
the analysis. There is a latitude -dependence of the stationary 
wave that was removed by parameterization. Inclusion of this 
additional y- structure would probably destroy the antisymmetry 
that leads to the meandering of the time -dependent solutions 
and that leads to an expectation of the vanishing of the phase 
oscillation of the wavenumber 2 component in the middle of the 
channel, a feature which is not observed in the experimental 
data. The results of LKP also suggest the possible presence of 
stationary wavenumber 4 and wavenumber 6 components. The solu- 
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tion of the time -independent equations should be extended to 
include these. 

It would be interesting to investigate the robustness of 
these conclusions in a nonlinear study. One question is 
whether one must specify the lower-layer velocity in the basic 
state. A nonlinear model is capable of generating a zonal flow 
correction. Nonlinear interactions could lead to the excita- 
tion of a topographic instability. It is desirable to see what 
sort of modes survive, their time evolution, and the types of 
final states which are present in various regions of parameter 
space . 



REFERENCES 


Blackmon, M.L., 1976: A climatological spectral study of the 

500 mb geopotential of the Northern Hemisphere. J . Atmos . 
Sci. . 33. 2285-2291. 

Blackmon, M.L., J.M. Wallace, N.-C. Lau, and S.L. Mullen, 1977: 
An observational study of the Northern Hemisphere Winter- 
time Circulation. J. Atmos. Sci.. 34. 1040-1053. 

Charney, J.G. and J.G. DeVore, 1979: Multiple flow equilibria 

in the atmosphere and blocking. J . Atmos . Sci . . 36 . 
1205-1216. 

Charney, J.G. and D.M. Straus, 1980: Form-drag instability, 
multiple equilibria, and propagating planetary waves in 
baroclinic, orographically forced planetary wave systems. 

J . Atmos . Sci . . 37 . 1157-1176. 

De Szoeke, R.A. , 1975: Some effects of bottom topography on 
baroclinic instability. J. Marine Res.. 33. 93-122. 

De Szoeke, R.A. , 1983: Baroclinic instability over wavy topog- 
raphy. J. Fluid Mech. . 130 279-298. 

Durney, B.R., 1977: The influence of mesoscale topography on 

the stability and growth rates of a two- layer model of the 
open ocean. Geo. Astronh. Fluid Dvn. . 9. 115-128. 

Fredericksen, J.S., 1979: The effects of long planetary waves 
on the regions of cyclogenesis: Linear theory. J . Atmos . 

Sci.. 36. 195-204. 

Gates, W.L., 1961: Static stability measures in the atmosphere. 
J. Met. . 18. 526-533. 

Hart, J.E., 1979: Barotropic quasi- geostrophic flow over aniso- 
tropic mountains. J . Atmos . Sci . . 36 . 1736-1746. 

Hart, J.E., 1981: Wavenumber selection in nonlinear baroclinic 
instability. J . Atmos . Sci . . 38 . 400-408. 

Li, Q.G., R. Kung, and R.L. Pfeffer, 1986: An experimental 

study of baroclinic flows with and without two -wave bottom 
topography. (In press). 


78 


79 


Nathan, T.R. , 1985: Interactions between unstable baroclinic 

and resonant, topographic waves in a two- layer model of the 
atmosphere. Ph.D. dissertation, State Univ. of N,Y. at 
Albany , 106pp . 

Niehaus, 1980: Instability of non-zonal baroclinic 

flows. J . Atmos . Sci . . 38 . 974-987. 

Pedlosky, J., 1979: Geophysical Fluid Dynamics . Springer- 

Verlag, 624 pp. 

Pedlosky, J., 1981: Resonant topographic waves in baro tropic 
and baroclinic flows. J . Atmos . Sci . . 38 . 2626-2641. 

Peixoto, J.P., B. Saltzman and S. Teweles, 1964: Harmonic ana- 

lysis of the topography along parallels of the earth. J . 
Geophvs . Res . , 69 . 1501-1505. 

Tibaldi,S. and A. Buzzi, 1983: Effect of orography on Mediter- 
ranean lee cyclogenesis and its relationship to European 
blocking. Tellus . 35A. 269-286 . 



APPENDIX A: 


ORTHOGONALITY CONDITION FOR ASYMPTOTIC MODEL 

I 

I 

The orthogonality condition will be presented for the case 
of U 2 * 0 and U;^ * U, but the result for general U 2 , to which | 

the method readily extends, will also be stated. The leading | 

I 

order problem from eq. (3.1) can be written in terms of the 

lower- layer streamfunction; i.e. in the following form: j 

i 

jr(0) - 0, (A.l) , 

with boundary conditions 

^ 2 ^®^ - ® at y - 0,1. (A. 2) i 

i 

i 

i 

The solution to this problem was given in eqs. (3.4)-(3.5), and < 

for U 2 “ 0 , it is given by 

- exp(ikx)sin Z-ny , (A. 3) 

— k/2{c+-ia' ± [c.^-2iac^ - a'^ + 4ot]^/^), (A. 4) 

where c^, c., o', and a were defined in eqs (3.6). 

The order hQ/F problem is given by representing eqs . 
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(3.10-3.11) in the form 

+ 8L]_(°)[sin2x - 0 (A. 5) 

with boundary conditions 

^> 2 ^) - + L 2 ^^^V> 2 ^^^ + 8L]_(°)[sin 2x V'2‘-°^y]-0 

(A. 6) 

In order to find , we apply the solvability condition 

by multiplying (A. 5) by the adjoint of (A. 3) and integrating 
over the domain: 

1 2jt 1 2;r 

/ / <^y+/ / (v>2(0)+)*j:(0)^2^^^ '^y 

0 0 0 0 

1 27T 

S 2x^2^^^y] 0. (A. 7) 

0 0 ^ 


where the superscript + represents the adjoint. 

The adjoint problem for ^ 2 ^^^ yields 

ir(0)+ - 0, (A. 8) 

where • V*2^^^ * exp(ikx) sin (i^ry) and is the same 

as except that the complex conjugate of , denoted 

c^(0)* is used in place of . 


The third term is zero since 
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L]^^°^[sin 2x » ( )exp(i(k+2)x) + ( )exp(i(k-2)x) 

j 

and therefore the product of this and exp(-ikx) integrates 
over X to zero. The first term in Eq. (A. 7) will be written 

1 

I 

/ / dx dy - (A. 10) I 

0 0 I 

The second term in Eq. (A. 7 ) is evaluated using a Lagrange , 

identity; i.e. , . 

0 0 0 0 I 

-/ [-iw(0)-U3x] [-iw^°^] [ ] I 

0 y-0 i 

(plus other terms which are zero when the b.c. are applied) i 

(A. 11) 

But the remaining boundary term integrates to zero when 

is expressed in terms of using the boundary i 

condition | 

4 

I 

- -4 sin 2x 

and substituted into the right-hand side of (A. 11). 

Also the second term on the left-hand side of (A. 11) is 
zero, leaving 


I - 0 


(A. 12) 
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In general, I is nonzero, implying - 0, This 
implies also that - 0, and then (3,10-3.11) become a 
forced differential equation for ^ 2 ^^^ terms of • 


” - 8 L 3 ^^®^[sin 2x (A. 13) 

i 

I 

I The solution ^ 2 ^^^ form 


^ 2 ^^) - g2+(y)exp(i(k+2)x) + g 2 _ (y)exp(i(k-2)x) (A. 14) 


The 0(ho/F)2 problem is written 


[sin2x ® (A. 15) 


with boundary conditions 

► 

i (A. 16) 

^ At 0(hQ/F)2 the orthogonality condition on (A. 15) yields 

I 

[' 

/ J dx dy +J / (^2^^^^)**^^^^ ” 

0 0 0 0 

i 1 2x 

J 2x^2^^V) (A. 17) 

^ 0 0 


yielding 


84 


w(2)i - I' + 47r2lF(w(®^-Uk+iri) [cos ixy (g2-y-g2+y)] |. 

y-0 
(A. 18) 

where the integral I is given by 
1 lit 

0,(2) I » J J (A. 19) 

0 0 

- ^ Mki 

where 

Mki ^ -(K2+F) [w(®)(K2+F)+2^k+(ri+r2)iK2] + F^ (2w((>) -Uk) 

(A. 20) 

The integral I' is given by 

I' - 47r2iF{a2+^(A++)-a2.?(A.+)+a4+|(A+.)-a4.e(A..) 

- l/ 2 (a 5 +-a 5 .) [^'+i2^2 (w(0) -Uk+iri) ] ) , 

(A. 21) 

where 

fi> s (w(0)-Uk)(k2+F) + k(/3+FU) + r^ik^ , 
and when i is even, 

^(A) a 2A[/i'+A2(w(0)-Uk+ir]^)](cos(A/2))/(i2^2.;^2)^ 

and when i is odd, 

?(A) a -2A[/i'+A2(«(0)-Uk+iri)] (sin(A/2))/(i2x2-A2) . 
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For general U 2 , equation (A. 18) becomes 

I - I’ + I" . (A. 22) 

where I ' and I ' ' are given by 

I' - 47r2i{(a)(*^)-Uik+iri)(b2.l2,4+'^2+^2,4-'*‘®2-^2,2+ 

-B 2 +I 2 , 2-+*^2-^^l , 2+'^2+^^l ,2-1 

- («(0) -U2k)F[bi. Ii ^ 2+'°l+^l , 2-+®l- ^1 , 0+*°l+^l , 0- 

+bi.Fl2^0-l-'^l-^^2,0-l 
- 4 ' tb2-l2,2+'’^2+l2.2-+®2-^2,0+*^2-i-^2,0- 

+ b2_FIi o+-b2+FIi,0-n> (A. 23) 


and 

I" - 4ff2i(„(0).xjj^k+r^i) . { 


d2 


d2 

b2- + B 2 . 

S2+y • 

b2+ + ^2+ 

dy2 

dy^ 


y-1 

+ b 2 .Fgi+y - b 2 +Fgi.y) cos i?ry | (A. 24) 

y-0 

and where 

1 

In,m± - f sin £iry — gn+(y) dy, 

0 dy^+i 

h± ~ -biT[(k±2)2-4]-Hb2q:F. 

B2± - -b2q:[(lct2)2-4]+(blq^^l)F. (A. 25) 



APPENDIX B: 


CALCULATION OF MERIDIONAL STRUCTURE COEFFICIENTS 

The functions g 2 + are given in eqs . (3.17) and are 

repeated here: 

S2+ ” ^1+ (A++(y-l/2)) + a 2 + sin (A++(y-l/2) + 

a 3 + cos(A+. (y-1/2) ) + sin( A+. (y- 1/2) + a54.cosi7ry 

(B.l) 

g2- “ ^1- cos (A.+(y-l/2)) + a 2 _ sin (A.+(y-l/2)) + 

a 3 _ cos(A..(y-l/2)) + sin( A. . (y- 1/2) ) + a^,cos£^y 

(B.2) 

First the coefficients a 5 + are found by the method of 
undetermined coefficients; i.e., we write a particular solution 
of eq. (3.14) in the form 

gp - ^ 5 + cos iiry, (B.3) 

resulting in 

as± — 03+ / + ci 2 ±] . 

(B.4) 
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These are substituted into the boundary conditions ( 3 . 15 ) 
to obtain, in the case of g 2 +. 

a]^^.cos(A^_^/2) - a2+sin(A4^./2) 

- a3^.cos(A+./2) - a4+sin(A^../2) + a54_ =■ 0, 

(B. 5 ) 

ai+cos(A^_^/2) + a2+sin(A^_j./2) 

- a3+cos(A+./2) + a4+sin(A+./2) + (-1) -a5+ - 0, 

(B. 6 ) 


ai+(-A^_4.+04)cos(A++/2) - a2+ (-A^_^+a4) sin(A^_^/ 2 ) 

- a3+ (-A+.+a4) cos(A+./ 2 ) - a4+ (-A+.+a4) sin(A+./ 2 ) 

+ (-i2^2+Q(^)a5+ - 05, 

(B. 7 ) 


a-l+ (-A+++04) cos(A++/ 2 ) + a2 + (-A+++ct 4) sin(A++/ 2 ) 

- a3^. (-A+.+a4) cos(A+./ 2 ) + a4+ (-A+.+a4) sin(A+./ 2 ) 

+ (-l)-«(-i27r2+a4)a5+ - (-D^aj. 

(B. 8 ) 

If Jl is odd, the sum (B. 5 )+(B. 6 ) and the sum (B. 7 )+(B. 8 ) 
yield only the trivial solution for a]^+ ^ 3 +- Then eqs. 

(B. 5 ) and (B. 7 ) are solved simultaneously to yield 
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&2+ sin(X++/2) - -X+.'^)+a^+, ' 

a 4 + (A++ 2 -A+. 2 ) sin(A+./2) - a 5 +(A ++2 .^ 2 ^ 2 ) .q,^^ 

If S. is even, the difference (B.5)-(B.6) and the difference 
(B.7)-(B.8) yield only the trivial solution for a 2 + and a 4 +. 

I 

Then eqs. (B.5) and (B.7) are solved simultaneously to yield 

! 

&l+ (A++2-A+.2) cos(A^_^/2) - a5+(A+.2.^2n-2).Q!5^, , 

a3+ (A++2.;^^_2) cos(A+./ 2) - a5+(i2rr2-A++2)+a5+. i 

i 

i 

j 

i 
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